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We present a general analysis of the role of initial correlations between the open system and an 
environment on quantum dynamics of the open system. 
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I. INTRODUCTION 

Proper understanding of quantum dynamics of open 
systems is a very important task in many areas of physics 
ranging from quantum optics, to quantum information 
processing and to quantum cosmology In general, 
one can assume an interaction between the open system 
denoted as A with the environment B. This environ- 
ment is a quantum system with the Hilbert space of an 
arbitrary dimension. The whole A + B system evolves 
unitarily. In most of the studies on dynamics of open 
systems it is assumed that the open system and its envi- 
ronment are at the initial moment of their joint evolution 
factorized |l],|2| , that is they are described by the density 
operator of the form 



Pab = PA® Pb 



(1.1) 



where pa is the initial state of system A and pb is the 
initial state of the environment. While the initial state 
of the open system A may vary the initial state of the en- 
vironment B is considered to be determined by external 
conditions. In this context it is natural to ask what is the 
time evolution of the open system Al Or in other words, 
what is the explicit form of the map %a ■ Pa — > Pa- I n 
order to answer this question one might follow the ar- 
guments presented in Ref. Q and to find the explicit 
expression for the density operator p' A : 

p A = Tr s {U A bPabU ab ) 

= J2^ Uab P a ® W ) U ab\h) 

= ^{p\^Uab\v)pa{v\^IJ\ b \\i) 

= Y, M ^PA M l», (I- 2 ) 



where 



M„ 



(\j\^VvU ab\v) 



(1.3) 



This is the well known Kraus representation j3j of a 
super-operator $^4 which has been studied and used in the 



literature broadly. On the other hand dynamics of open 
system in which initial correlations between the system 
and the environment is taken into account has not been 
analyzed in detail yet. Taking into account recent inter- 
est in quantum entanglement within the context of quan- 
tum information processing it is appropriate to study in 
detail the role of (quantum) correlations on dynamics of 
open quantum systems. Some particular aspects of this 
problem have been discussed earlier in Refs. . In par- 
ticular, in Ref. Q decoherence and initial correlations in 
quantum Brownian motion have been analyzed, while in 
Refs. H a motion of a free particle coupled to a lineary 
dissipative environment has been studied. The conection 
between the initial correlations and the complete positiv- 
ity has been discussed in Ref. P|. 

In the present paper we present a general analysis of 
the role of initial correlations between the open system 
and an environment on quantum dynamics of the open 
system. 



II. THE ROLE OF INITIAL CORRELATIONS 

In this section we will investigate the evolution of an 
open system A which is initially correlated with the en- 
vironment B. Let us denote by (7j the generators of the 
group SU(N) |t]] where N is the dimension of the Hilbert 
space of the system A. In addition we denote by Tj the 
generators of the group SU (M) where M is the dimen- 
sion of the Hilbert space of the environment B. Using 
this notation the most general density matrix of the sys- 
tem A and the environment B reads as 



1 



PAB 



NM 



(Iab + ® 1 B + 0jt; 



+ ~{ij<Ji®Tj) (2.1) 



while the density operator of the open system A is ob- 
tained via "tracing" over the environment 



Pa = Tfs (Pab) 



1 

N 



(1 A + OtiVi). 



(2.2) 



So let us assume that the state (2.1) is the initial state of 
the whole A + B system which evolves according to the 
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given unitary matrix Uab- Can we describe in this case 
the evo lutio n of the subsystem A in the form analogous 
to Eq. (1.2)? In order to answer the question we have to 
insert into the equation 



p' A = Tr s (Uab pab U\ b ) 



(2.3) 



the expression ( |2.l| ) for the density operator pab which 
results in 



Pa 



Y,{lA u AB jj{U + am) ® h{l B + PjT 3 )u\ B \p) 

J2(p\U AB {pa®Pb) U A b\p) 
+ {llJ N T ]) ^ ® T i U ab\p) (2-4) 



7VM 

After rewriting 7^ as 7y = N.M 7^ 
from Eq. (|2.4|) the expression 



Qti/3j we obtain 
p' A = Y J M tiV p A Ml v (2.5) 

7-j-CTj (g) (Tj 



where the operators -M^ are given by Eq. (1.3). We see 
that the resulting density operator describing the open 
system A during the time evolution consists of two terms. 
The first term corresponds to the standard Kraus rep- 
resentation with no initial correlations as discusse d in 
Sec. |. The second term in the right-hand side of Eq. ( |2~5| ) 
depends only on the correlation parameters 7^ which do 
not explicitly depend on the particular choice of the initial 
state of the open system A (see below). In other words, 
these parameters cannot be determined by performing a 
local measurement on the initial state of the system A. 

This second term makes the equation (2.5) inhomoge- 
neous though linear - we will discuss this consequence of 
initial correlations between the system and the environ- 
ment in the following section. 

Example 1. We may regard 7^ ■ as additional parame- 
ters which together with the initial state of the environ- 
ment ps and the unitary operator Uab determine the 
time evolution of the open system initially prepared in 
the state pa- 

To illustrate the possible significance of 7^ we will 
study a simple model describing dynamics of two qubits 
(spin- 1/2 particles). In this model one of the qubits (A) 
plays the role of the open system while the second qubit 
(B) plays the role of the environment. Let the unitary 
evolution operator U A b acting on the joint system of 
these two qubits is given by the expression 



where H is the Hamiltonian 



H = a x <g> -(1 - <y z ) + 1 <g> hi + g z ) 



(2.7) 



with <jj being Pauli matrices. The interaction described 
by the Hamiltonian (12/71) corresponds to the well-known 
controlled-NOT gate gj. 

Let us consider two initial conditions p AB and p AB for 
the two-qubit state, which in the computer basis {|0), |1)} 
read 



o {1) 
Pab 

S 2 



! |00)(00| + |/3| 2 |11)(11|, 



Pab = (a|00) + /3|ll»(a*<00| 



(2.8) 



Obviously, the qubits A and B in these two-qubit states 
are in the same state, i.e. 





= Tr B 




= Tr B 


(2) " 

Pab 




= Tr A 


'& 


= Tr A 


(2) " 

Pab 



(2) 

Pa 

(2) 

Pb ■ 



(2.9) 



On the other hand the parameters 7^ in the two-qubit 
states (2^8) are different. Therefore, this simple model 
with identical states of subsystems but different correla- 
tions will illuminate the role of the correlations on dy- 
namics of open quantum systems^ 

With the unitary evolution @.6| ) the two-qubit sys- 
tems with the two initial conditions (2.8) evolve at time 
t = it/2 into states such that the system A is described 
by the two density operators 



^(t = 7r/2) = i(l + ff3 ), 



1 



^(t = 7T/2) = -[l + (M S 



>3] 



(2.10) 



U = e- lHt = Icost 



iH sin t . 



(2.6) 



respectively. We stress here that the open system has 
been in both cases in the same initial state, i.e. p^ — 

p { 2 ] = H 2 |0)(0| + |/3| 2 |1)(1|, the environment itself was 
in both cases initially in the same state as well. But due 
to different initial correlations between the system and 
the environment the open system has evolved into two 
different states p A 1 \n/2) and p A \n/2). 

This example illustrates that the initial correlations 
between the system and its environment may play impor- 
tant role in the dynamics of open systems. Moreover, in 
most of physical situations such correlations are present 
and therefore they have to be taken into account. 



III. MASTER EQUATION 

As follows from our previous discussion both the state 
of the environment and the initial correlations between 
the environment and the open system play significant 
role in the dynamics of the open system. Therefore, in 
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order to characterize completely the evolution, it is nec- 
essary to determine (fix) the set of the parameters 
i.e. the state pb of the environment, and the parameters 
{7-} describing the correlations. The parameters o>i, (3j 
and 7L are arbitrary conditioned that the matrix pab 
describe a real physical state of the system AB, that is, 
it is a density matrix. Specifically, if we represent one 
particular choice of parameters {o^, (3j, 7-} as a point in 

a (N 2 ■ M 2 - 1) dimensional space R(^ 2 -* f2 -i) ) then the 
set of physically relevant parameters {o^, 7-} form a 
convex subset S in the space J\( n2 -M 2 -i) _ For example, 
in the case of on ( the same holds for f3j and 7^ ) there 

is only a subset Oa in the space from which we 

can chose the parameters {cti} so that the density matrix 
Pab is a reasonable density matrix. Moreover this subset 
Oa depends on the choice of the remaining parameters 
{Pj} and {7-}. In other words the subset Oa is different 
for different choices of {/3,} and For example, if 

the system A is in a pure state then the only possible ini- 
tial density matrix of the system A and the environment 
must have the form pab = PA® Pb, so that all 7^ have 
to be zero. Or equivalently if some of the parameters 7^ • 
are not zero then the state pa can't be a pure state. 

Sometimes it's very useful to describe the evolution of 
the open system in a form of a master eq uati on. In order 
to do so we firstly rewrite the evolution (2.5) in terms of 
the left-right super-operator acting on the density oper- 
ator pa 



p A (t)=T{t) p A (0) + £(t), 



(3-1) 



where £(t) is the inhomogeneous term which has its origin 
in the presence of initial correlations between the open 
system and the environment, i.e. from Eq. (2.5) we have 



£(*) = !>! Cab 



7»- 



7 3 u ab\v)- 



A" J 



We stress once again that the operator £(t) does not de- 
pend explicitly on the initial state of the open system A, 
only the range of possible values of correlations is deter- 
mined by the choice of pa and ps (see the discussion 
above). As follows from Eq. ( |2.5| ) the left-right action of 
the super-operator T{t) is equal to the following normal 
action 

f(t) pa(o) = J2^pa(o)mI v . 



^From our previous comments it follows that the choice 
of the initial correlations restricts a set of density oper- 
ators pa for which Eq. ( |3.l[ ) can be used. For instance, 
for pure states the term £(t) is always zero. Therefore, if 
we would use Eq. ([O]) with nonzero £(t) for describing 
dynamics of an open system initially prepared in a pure 
state, we might end up with a completely unphysical sit- 
uation. As discussed above this subset is determined by 



the condition, that dynamics (1.2) has a physical mean- 
ing. This restriction reflects quantum nature of correla- 
tions between the system and the environment and have 
to be taken into account in the derivation of dynamics 
of open quantum systems which are initially correlated 
with the environment. 

We have to keep in mind that there is always only a 
subset Oa of all de nsit y matrixes of the system A for 
which the equation (3.1) with a given £(t) is valid. If, for 
example, £(t) = then the equation (3.1) is valid for all 
Pa and Oa — Sa where Sa is a set of all density ma- 
trixes of the system A. Unless £(i) equals to zero Oa is 
a subset of Sa- 

After this preliminary comments we derive the master 
equation following the formalism presented in Ref. [||. 
Differentiating Eq. (3.1) according to time we obtain 



|paw = |f ( t) ,40) + 1 m 



(3.2) 



When we substitute pa(0) which formall y ca n be deter- 
mined with the help of Eq. (|1]) § into |T|) we find 

d . . / d ~, . 



J- [ PA (t)-m + §- t m (3.3) 



If we introduce a notation for the Liouvillian super- 
operator 



X 



d it., \ 1 



(3-4) 



then the master equation can be rewritten in the follow- 
ing form 



§- t -X) [ /M , /)-;(/)]=,» 



(3.5) 



If the initial correlations were zero, then the master 
equation ( |3.5| ) reduces to the well known form (see for 
instance, Ref. ||) 



(3.6) 



where the operator X is the same as in Eq. ([T^). Taking 
into account the fact, that £(i) does not depend on the 
initial state pa{0) we can introduce the operator 



(3.7) 



and rewrite the master equation (3.5) in an inhomoge- 
neous form 



d_ 

dt 



X\ p A (t)=T{t), 



(3.8) 



We stress once again that the super-operator X depends 
only on the initial state of the environment Pb and the 
parameters of the unitary evolution Uab , while the whole 
information about the initial correlations between the 
open system and the environment is in the operator ^(t). 
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IV. DISCUSSION 

Till now we have studied how initial correlations be- 
tween the open system and the environment can influence 
the time evolution of the open system. We have found 
that these correlations play an important role which can- 
not be neglected. In this section we will investigate prop- 
erties of super-operators (evolutions) $a acting on an 
open system which is a part of the composite system 
(open system + environment). It is assumed that two 
parts of the composite system can be initially correlated. 
The composite system is considered to be closed so that 
it evolves unitarily according to a given unitary operator 
Uab- In what follows we will assume this evolution of 
the "Universe" to be given. 

Firstly we define the most general super-operator (evo- 
lution) which originates from a given Uab- 

Definition IV. 1 Let A is the system of interest, B is 
the rest of the Universe (the environment), and Uab * s 
a given unitary evolution on the whole system. Let us 
consider a map 



V ■ Pa~> Pab 



(4.1) 



which means that for each pa we choose one pab from a 
set of all possible density matrices of the Universe such 
that 



Tib (pab ) = Pa 



(4.2) 



The super-operator which describes the most general 
evolution of the system A is given by the expression 



PA ->■ p A 



Pa = Tr B (u A bPabU ab 



(4.3) 



(4.4) 



The map V in Eq.(f.l) is related to the preparation of 
the state pa of the system A. We note that while prepar- 
ing the state pa of the system A the state of the Universe 
is changed as well. That is, in any act of the preparation 
of the system A we prepare a state pab rather then an 
isolated state pa of only the system A without affect- 
ing the system B. For this reason pab describes a (cor- 
related) state of the open system and the environment 
(Universe). Moreover, since the preparation is an act in 
which a classical information is encoded into a quantum 
system the map V is not necessarily linear. Therefore the 
state pb = Tya[pab\ might depend (even in a nonlinear 
way) on the state pa- For instance we can imagine the 
map V of the form: V(pa) — Pa® Pa which describe the 
action similar to quantum cloning which obviously is not 
possible within the framework of linear quantum mechan- 
ics. But can easily be performed at the level of prepara- 
tion of quantum states. Analogously we can imagine a 
map V(pa) = Pa® Pa, where p T is a transposed state. 



Taking into account that Uab is fixed then the only "free- 
dom" in controlling the dynamics is the choice of the map 
V. 

It is clear from the construction that the super- 
operator $ is a trace preserving map and that the final 
operator p A is Hermitian and positive, i.e., it is a valid 
density matrix. In what follows we will study some as- 



pects of the of the evolutions of the form [IV. 1 
a 



L 1 



From the definition IV. 1 it follows that for a 



given Uab and an arbitrary map V not all evolutions $^4 
can be realized. On the contrary there exists Uab and V 
such that a given $^4 can be realized. To see this let us 
consider a following example. 

Example 2. Using the scenario (4.4) we can perform 
any map $ : pa — > p\ on a given (known) initial state 
PA of the system A. Specifically, let $ : pa — > p'a is a 
given map. We assume that the map V acting during 
the preparation of the system A is such that the compos- 
ite system has been prepared in the state pab 

Pab — PA® Pb 

such that pb = p'a The unitary transformation which 
realizes the desired map is then taken to be 



Uab =J2\ i )A{j\®\j) 



Obviously, there is nothing surprising here since if we 
know the initial state of the system pa exactly then we 
can perform an arbitrary map on the system. In some 
sense this situation corresponds to a classical physics 
when a complete knowledge about the state of the system 
is implicitly always assumed. Knowing the initial state 
precisely we can perform any map we want. 

b. Till know we had not considered the linearity con- 
dition in association with the evolution As we have 
already commented the unitary evolution Uab is by the 
definition linear, but the preparation map V might be 
nonlinear. At this moment we can ask which conditions 
on V has to be imposed so that $^ is linear. In order 
to proceed we remind us the definition of the linearity of 
the evolution $^4. Namely, $,4 is linear if 



Now it is clear that if V is linear, in a sense that 



(4.5) 



V 



A I 



(4.6) 



then the evolution $^ is linear. The linearity of V is 
a sufficient condition for the linearity of $,4. On the 
other hand it is not the necessary condition. We might 
imagine a nonlinear map V such that %a is linear. To 
understand this we formally represent the evolution $^4 
as %a — TybUab, where we use notation such that 
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Uab(pab) = UabPabU\ b . 
(4.5) can be expressed as 



Then the linearity of 



Tr bUabV 




A = Y, X ^aUabV{ P 



(0\ 
A ) 



(4.7) 



= Tt b Uab 



^From this last equation it follows that if the map V 
is linear, then %a is linear as well. On the other hand 
from the linearity of $,4 does not follow that V is linear. 
This is a consequence of the property of the part ial trace 
operation Tre. Specifically, from the identity (4/7) the 
equality 




(4.8) 



does not follow. 

c. Next we will consider consequences of another 
possible restriction on V . Namely, let us consider a rather 
frequent condition, that the state of the environment ps 
does not depend on the state of open system pa- That is 



Tr A V(pA) = Pb = const 



(4.9) 



MP 



for all pa- If Pa is pure, then under the condition 
the map V is uniquely defined such that V(pa) = 
On the oth er hand if pa is impure, then under the con- 
dition (19) the map V might not be uniquely specified, 
i.e. correlation bet wee n A and B can play a role. 

If the condition (4.9) is fulfilled and in addition the evo- 
lution $a in the definition IV. 1 is linear, then the map 
V can be chosen such that V(pa) — Pa® Pb for all pa- 
But this means that the evolution %a can be represented 
in the Kraus representation |l4| ] . Consequently, this map 
is completely positive M. 
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